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We propose a quantum pump mechanism based on the particular properties of graphene, namely 
chirality and bipolarity. The underlying physics is the excitation of evanescent modes entering a 
potential barrier from one lead, while those from the other lead do not reach the driving region. 
This induces a large nonequilibrium current with electrons stemming from a broad range of ener- 
gies, in contrast to the narrow resonances that govern the corresponding effect in semiconductor 
heterostructures. Moreover, the pump mechanism in graphene turns out to be robust, with a sim- 
ple parameter dependence, which is beneficial for applications. Numerical results from a Floquet 
scattering formalism are complemented with analytical solutions for small to moderate driving. 
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I. INTRODUCTION 

Ratchets and pumps are devices in which spatio- 
temporal symmetry breaking turns an ac force without 
net bias into directed motionP^ If the time-dependence 
enters via only one parameter, the pump current van- 
ishes in the adiabatic limit. 3 Therefore, single-parameter 
pumping requires non-equilibrium conditions enforced by 
driving beyond the adiabatic limit. This distinguishes 
single-parameter pumps from devices that operate with 
oriented work cycles, like turnstiles^ or sluicesP Quantum 
pumps can be implemented with quantum dots in a two- 
dimensional electron gas (2DEG) driven by microwaves^ 
surface- acoustic waves P ac gate voltaged^^ or non- 
equilibrium noiseP^ Here we propose a single-parameter 
pump based on the particular properties of graphene and 
show that these display a broad-band response, in con- 
trast to 2DEG-based ratchet s or pu mps for which isolated 
resonances govern the effect JsE^EsD 

Owing to the chiral and gapless nature of charge car- 
riers in graphene, an electron hitting a potential step 
in this material may prop agate forwards as a hole with 
opposite momentum! 14 * 15 * Due to this Klein tunneling it 
is difficult to confine electrons by electrostatic potentials. 
This phenomenon even extends to evanescent modes, i.e., 
modes that decay exponentially as a function of the bar- 
rier penetration. In graphene, electrons populating such 
modes can tunnel a large distance. 16 This impediment 
to electrostatic confinement constitutes a drawback for 
many switching and sensing applications. 

An inspection of present quantum pump designs sug- 
gests that graphene pumps would be negatively affected 
by this issue as well. For example, pumps have been real- 
ized with e lectr ostatically defined double quantum dots 
in a 2DEGP^ In these experiments, a driving field in- 
duces dipole transitions of an electron from a metastable 
state below the Fermi energy in the, say, left dot to a 
metastable state in the right dot. Subsequently the elec- 
tron will leave to the right lead, and an electron from 
the left lead will fill the empty state in the left dot. The 



emerging pump current therefore requires resonant con- 
ditions with a pair of energetically well-defined states, 
i.e., well isolated states with long life times, which do 
not exist in graphene. 

Nevertheless, the alternative, graphene-specific mech- 
anism identified here allows realizing highly efficient 
pumping, which may indeed outperform conventional de- 
vices. The large pump current emerges in the bipolar 
regime around the Dirac point. This is due to a scattering 
process where a whole continuum of evanescent modes is 
promoted into unidirectionally propagating states, which 
couple well to the leads because of chirality. Since this 
mechanism does not rely on intricate resonance condi- 
tions, but invokes graphene's intrinsic features of bipolar- 




FIG. 1. Sketch of a single-parameter graphene pump con- 
nected to two leads. The left half of the graphene ribbon is 
exposed to an ac gate voltage. Dispersion of (b) graphene 
and (c) a 2DEG, with longitudinal and transverse momenta 
k and q, including the branches of evanescent modes. Photon 
absorption (and in the case of graphene also emission) may 
excite an evanescent mode to propagating. 
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ity and chirality, the pump current is robust and displays 
a simple parameter dependence curve, as is desirable for 
electronic sensor applications. 

This paper is organized as follows: In the next section 
we describe our proposal for a single-parameter pump 
and give typical values for length and energies scales. 
In Sec. |III| we introduce the Floquet scattering theory 
that applies to non-equilibrium pumping in graphene. 
We first present the general formalism and then con- 
sider the weak driving regime, where we derive analytical 
expressions for the one-photon transmission probability. 



In Sec. IV we present numerical results for the pumped 
current comparing the performances of graphene and a 
2DEG pump. Moreover, we compare the numerics with 
the semiclassical approximation derived in the preceding 
section and find good agreement. The specific pump- 
ing mechanism by excitation of evanescent modes is ex- 
plained in Sec. |IV A| Finally, we conclude in Sec. |V| 



II. DESCRIPTION OF THE SYSTEM 

As a single-parameter pump setup, sketched in Fig. [T^l, 
we consider a graphene ribbon of length 2L and width W 
attached to two metallic electron reservoirs, with its left 
half driven by a time-dependent gate. The driving by 
variation of only one parameter excludes the emergence 
of an adiabatic pump current However, considering 
that spatial symmetry is broken by the placement of the 
gate, a finite dc current arises in non-equilibrium condi- 
tions, which can be achieved by non-adiabatic driving. 
In order to assess the importance of graphene's chirality 
and bipolarity in the non-adiabatic context, we contrast 
our results with those found for the corresponding setup 
of a 2DEG in a semiconductor heterostructure. 

We assume that the sample sizes are smaller than the 
mean free path and that W ^> L, so that the effects 
coming from the boundaries, electron-electron interac- 
tions and disorder play a minor role. Besides, we re- 
strict ourselves to the low energy regime such that the 
Dirac approximation remains valid. Then, our quasi one- 
dimensional model contains three natural energy scales, 
namely the driving frequency lj, the driving amplitude t/, 
and the energy associated with the length L of the de- 
vice. Additionally, the leads' Fermi momentum hk^ ^ 
becomes a relevant scale in a 2DEG, since the con- 
tact resistance depends on its value. For graphene with 
Fermi velocity vp ~ 10 6 m/s and L = 5/im, the lat- 
ter is E^ = hvp/L « 0.13meV. For a 2DEG with 
the same geometry and effective mass m* = 0.067ra e 
(GaAs/AlGaAs), it is roughly four orders of magnitude 
smaller, = h 2 /2m*L 2 « 0.02/ieV. Driving beyond 
the adiabatic limit requires frequencies uo > El/H. Elec- 
trons can then absorb or emit photons, and after a short 
transient period, these excitations will establish a non- 
equilibrium population of the electronic states. 



III. FLOQUET SCATTERING THEORY IN 
GRAPHENE 

A quantitative description of nonadiabatic charge 
transport is provided by Floquet scattering theory. Sub- 
sequently, we show that this theory provides the proba- 
bility Tj^i for an electron to be scattered from the left 
to the right lead under the absorption or emission of n 
photons, and with it we can calculate the dc current.^ In 
order to focus on the graphene-specific features, we will 
restrict ourselves to the zero temperature limit so that all 
electronic states below the Fermi energy Ep are initially 
occupied. 



A. General formalism 

When scattered at a periodically time-dependent po- 
tential, an electron with initial energy e may absorb or 
emit \n\ quanta of the driving field (n < corresponds to 
emission), such that its final energy is e + nftw. This is 
embodied in the decomposition of the transmission prob- 
ability from the left to the right reservoir 



T LR {e)= £ Tj$(e) 



(1) 



The corresponding time-average d dc c urrent is given by 
the generalized Landauer formuld^^l 



ge 
h 



J deJ2 [Tii(e)h(e)-T$(e)f R (e)], 



(2) 



where e is the electron charge and /(e) is the Fermi-Dirac 
distribution. Spin and valley degeneracy of graphene 
is responsible for the pref actor g = 4, while g = 2 
accounts for the spin in a 2DEG. The application of 
Eq. pi to graphene (or any other two-dimensional ma- 
terial) requires extending the summation to transverse 
momenta hq. For a driving field that breaks reflec- 
tion symmetry, one generally finds Trl(c) ^ Tlr(c). 
Then, even when the leads are in equilibrium such that 
/z,(e) = = f( e )i a ne t current may flow, and a 

pump current emerges, 



= de/(c)AT( C ), 



(3) 



where AT = T LR - T RL . 

For the computation of the transmission probabilities 
we adopt the Floquet scattering formalism of Ref.l20land 
consider electrons in two dimensions under the influence 
of a time-dependent potential described by the Hamilto- 



nian 



H(t) = H (x) + U(x) cos(cjt). 



(4) 



Here, Hq comprises the kinetic energy and the static 
potential V(x), while U(x) is the profile of the time- 
dependent potential with frequency uo. Since the po- 
tential is ^-independent, the transverse momentum is 
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FIG. 2. Quasi-one-dimensional scattering potential. Profile 
of the pump, modeled as a barrier for which the left half 
(region I) experiences an oscillatory gate voltage. Regions L 
and R correspond to highly doped leads. Dotted lines denote 
sidebands to which the electron energy changes by absorption 
and emission of photons from the driving field. 



conserved and the problem becomes effectively one- 
dimensional. We assume that both the static potential 
V(x) and the driving profile U(x) are piecewise constant, 
and that U(x) = outside the scattering region; see 
Fig. [2] For graphene, Hq = hvpk • a + V(x) with the 
wavevector k = ±ke x + qe y . The free solution with en- 
ergy E reads 



by a separation ansatz. The resulting solutions 

4>± e (x,t) = e- m ^ u ^ si ^^f +nhu _ Ve (x) (6) 

CO 

= ^^(Pf/Me-^VU-^W. 



n — — 00 



comply with the requirement of time-periodicity provided 
that the separation parameter n is of integer value. The 
separation parameter labels all possible solutions and de- 
termines the time- averaged energy e + nhuj. The Bessel 
function of the first kind J n stems from the relation 
exp[— iz sin(co>t)] = ^2 n J n (z) exp(— inuot). Note that in 
the present case, XJi is non-zero only in the driving re- 
gion t — l. The ansatz (|6| has also been used to study 
photo-assisted tunneling in graphene! 21 * 22 * 

Condition (ii) means that in region L, the Floquet state 
consists of an incoming plane wave and a reflected part, 
(j) L = vt-v L +X) n r n^n,L' wnile in region R, we have only 
an outgoing state, <j)R = ^n^^n r- With the normal- 
ization chosen for the free solutions ([5| , the coefficients of 
the latter superposition relate to the left-to-right trans- 
mission probability under absorption or emission of n 
quanta according to 



T^(e) = \t n (e)f 



(7) 



y/2\E\k/K 



\E\/hv F \ 

±k + iqj 



(5) 



where k is positive and fulfills the dispersion relation 
E 2 = H 2 v 2 7 (k 2 + g 2 ), while ± is the sign of the corre- 
sponding current. The normalization has been chosen 
such that propagating waves have unit longitudinal cur- 
rent, vpi^y (J x ( Pe = il- This is convenient since then 
the coefficients of a superposition become probability am- 
plitudes. For evanescent solutions with imaginary longi- 
tudinal wavenumber k = in, the current vanishes. 

According to the Floquet theorem, the Schrodinger 
equation with a time-periodic Hamiltonian H(t) = H(t + 
2tt/uj) possesses a complete set of solutions with struc- 
ture ip = e~ %et l h (\)(t\ where the Floquet state <j>(t) = 
(/)(t-\-2ir/uj) obeys the time-periodicity of the Hamiltonian 
and e is the quasienergy. Here we are looking for Floquet 
scattering states, i.e., solutions of the Schrodinger equa- 
tion that (i) are of Floquet structure and (ii) have an 
incoming plane wave as boundary condition. For clarity, 
we derive here only the transmission from left to right — 
for the opposite direction, it follows by simple re-labeling. 

Condition (i) is equivalent to employing for the wave- 
function in any of the four regions £ = L^l^r^R the 
ansatz ipi(x,t) = e~ let / h <j)i(x,t). The time-periodic parts 
4>i(x, t) still have to be determined, while the quasienergy 
e turns out to equal the energy of the incoming wave. In- 
serting ip£ into the Schrodinger equation of region £ yields 
for (j>i(x,t) a partial differential equation which we solve 



In the scattering regions I and r, the Floquet solution 
must be a superposition of the states ([6|, i.e., (f)£ = 

Z) n ^Vn,€ + En r «Vn,£' Tne coefficients t$ and 
follow from the requirement that the scattering states 
have to be continuous at any time. These matching con- 
ditions can be written as a set of linear equations, with 
the incoming wave appearing as inhomogenity. In matrix 
notation it reads 



M , • (r / r (0 f (r) r (r) t A 



(8) 



where </5j in = pJ-Vl denotes the incoming wave written 
as a 6-dimensional vector. The 6x6 matrices M n ^ n > are 
constructed from the wavefunctions ([6| evaluated at the 
interfaces. They can be written efficiently as 



M(e - n'tuJ) Jw-n (Ui/hu) P £ , (9) 



where M(e) contains the wave matching condition at en- 
ergy e in the absence of driving. The matrix Pg is a 
projector to region constructed such that MPg con- 
tains only wavefunctions from region L In particular, 
Pi = diag(0, 1, 1, 0, 0, 0) projects onto the driving region. 
Equation ([8| fully determines the transmission and reflec- 
tion amplitudes t n and r n . The tight-binding version^ of 
this method is suitable for studying ac driving of smaller 
carbon-based conductors, such as ribbons with a size of 
only a few lattice constants^ or nanotubes ! 24 * 25 * Via time- 
dependent density functional theory, it can be generalized 
to the presence of interactions.^ 
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For electrons in a 2DEG with effective mass m*, the 
static Hamiltonian reads Ho = h 2 k 2 /2m* + V(x). Its free 
solutions are scalar fields, but since also their derivatives 
must be continuous, it is convenient to write them in 
spinor notation, 



<Pe 



±ik 



(10) 



with the dispersion relation E = h 2 k 2 /2m* . The normal- 
ization is such that the first vector component has unit 
current, ^*(<p% ^ (—ihd x )(pg ± = ±1. With these ingre- 
dients, the Floquet scattering formalism can be directly 
applied to the 2DEG case. 



B. Weak driving limit 

The set of linear equations ^ can be solved analyt- 
ically in the limit of small driving amplitudes or large 
frequencies, such that p = (U/2hu) 2 <C 1. We aim at 
finding analytical expressions for the one-photon trans- 
mission probabilities T^ 1 ) in terms of the static trans- 
mission for Klein tunneling at energy e with respect to 
the top of a high barriei^ of length L, 



T(e,q) 



fc e 2 + <Z 2 sin 2 (fc £ L)' 



(11) 



where k e = [(e/hvp) 2 — q 2 } 1 ^ 2 . The semiclassical limit 
L — )> oo of the above relation is obtained by inte- 
gration over fast oscillations, which yields T(e, q) « 
y/l — (hvpq/e) 2 . In order to relate these expressions to 
our formalism, we extract from the solution of Eq. (JsJ) in 
the undriven limit, U = 0, the transmission amplitude 
to by multiplication with the vector p^ = (0, 0, 0, 0, 0, 1) 
and obtain 



T{e) 



\^ R M-\e)vh 



(12) 



Next we simplify Eq. (|8| using for the Bessel func- 
tions the approximations Jq(U/Huj) = 1, J±i(U/Hlj) = 
±U/2Huj, while J n (U/Huj) = (U/Huj) for U/Huj <C 1. 
Thus, to first order mU/fiw, only the Floquet indices n = 
0, ±1 remain. The resulting 18 x 18 matrix (M n>n /) can 
be inverted via the approximation [A + (U/2ftw)B]~ 1 = 
A~ x — (U/2ftw)A~ 1 BA~ 1 , which provides t±\ and, thus, 



r (±i) 



^ R M- 1 (e T nw)M(e)P l M- 1 (e)^l h 



(13) 



The calculation of the elastic transmission T^ up to 
leading order in p is more tedious, but fortunately not 
required, because our Hamiltonian obeys time-reversal 
symmetry. Therefore, T^(e) = T^(e) which implies 
that the elastic channel does not contribute to the pump 
current.^ 

Our numerical calculations will reveal that the relevant 
contributions to the pump current stem from modes that, 



in the absence of driving, are evanescent in the barrier 
region. For such modes, Eq. (|13|) can be connected to 



Eq. (12) in a simple way, because the modified source 



term in Eq. (13) has an invariant forward amplitude, 
M(e)P l M~ 1 (eVtin = $~ >in + «"^Z,in- Although or is 
not necessarily zero, it does not affect the transmission, 
since it merely redefines the reflection amplitude r$. As 
a result, the one-photon transmission, Eq. (13), becomes, 



besides a prefactor p, identical to the static transmission 



( 12 ) in the evanescent region, 



T i L f ) (e)=pT(e±tuv,q). 



(14) 



For the evanescent modes entering from the right, the 
same reasoning lets us conclude that the waves decay 
exponentially and do not reach the driven region I. Con- 
sequently, PjM(~ 1 )(e)^ jin « 0, and T^ 1} (e) <C p can 
be neglected. Then the net transmission appearing in 
Eq. (|3| becomes 



AT = p [T(e + q) + T(e - fiw, q)} + 0(e 



-2kL 



). (15) 



The magnitude of the correction reflects the fact that 
we have neglected the exponentally small transmission of 
the evanescent modes which decays with the imaginary 
wavenumber ik = i[q 2 — (e / hv p) 2 } 1 ^ 2 • Since < T < 1, it 
follows that the maximal net transmission for evanescent 
waves is 2p. 



IV. NON-ADIABATIC PUMP CURRENT 

In the relevant weak driving regime U <C Huj, the 
absorption or emission probabilities are of the order 
p = (U/2hu) 2 < 1 [see Eq. ([HJ], so that AT - p. More- 
over, for short and wide systems (width W ^> L), the 
total current takes the form of an integral over modes 
q. Hence, the current ^ at a given Fermi energy Ep 
(measured from the Dirac point of the barrier) can be 
expressed as 



h \2huj J 



/E p pC 
de 
-oo J — c 



dq 



AT 



(16) 



We compute AT numerically by wave-matching in Flo- 
quet space [see Eq. Using typical parameters L = 
5/im, W/L = 4, U = 40/ieV, Huj 2meV (around 
500GHz), we obtain the results shown in Fig. [3] As 
long as p <C 1, the structure of the transmission for 
graphene depends only on the product Tcj, unlike for the 
2DEG. Graphene develops a far larger pump current than 
a 2DEG, saturating to an L- and cj-independent maximal 
value / max « 3nA for \Ep\ > fi^, while the smaller pump 
current in the 2DEG case already saturates at \Ep\ > 0. 
The semiclassical approximation described in Sec. |IIIB[ 
valid for huj ^> El, yields 



\e f \ 

Huj 



)f£, \E F \ <hw 



\E F \ > hoj 



(17) 
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2DEG 



FIG. 3. Pump current / for graphene and 2DEG. Dashed 
lines mark the semiclassical result, which closely matches the 
numerical ones. The upper inset is a blowup of the 2DEG 
case. The lower inset depicts the differential current demon- 
strating that the main contribution to the pump current arises 
from an energy range ftw around the Dirac point, populated 
by modes that are evanescent in the barrier region. 



where the energy E\y is the analogue of El with L re- 
placed by the width W. The corresponding calculation 
for a 2DEG provides the result 



In 



{U/2f 



^2k F ° o) WE^ 



0, 
1, 



E F < -hw 
-huo < E F < 
E F > 



These approximations are plotted as dashed curves in 
Fig. [3] Unlike for graphene, the pump current in the 
2DEG depends on the leads' large Fermi momentum k^^ 
(around 12 nm -1 for gold electrodes), which has been as- 
sumed equal for both leads. The relative pump perfor- 



mance v 



max 
G 



max 
N 



hkp°^ /m*VF-, assuming gold 



electrodes and GaAs/AlGaAs 2DEGs, is v « 20, i.e., the 
pump current in the graphene device is much larger than 
in the 2DEG device. 



A. Direction-dependent excitation of evanescent 
modes 

In the following we show that the large and robust 
pump current of the graphene device stems from a mech- 
anism whereby the ac field promotes evanescent modes 
from the left lead with probability p into propagating 
modes, unlike evanescent modes from the right lead, 
which couple poorly to the driving region. In order to 
support this picture, we consider the differential response 
dl/dEp, shown in the lower inset of Fig. [3] It indicates 
that the main contribution to the current stems from the 
bipolar regime |e| < hou in the case of graphene, or from 
the gap boundary region — thus < e < in the case of the 
2DEG. In the absence of driving, these energy ranges are 
populated by electronic modes that become evanescent 




FIG. 4. Scaled net transmission AT/p, where p — (U / '2hw) 2 , 
for (a) the graphene and (b) the 2DEG pump as a func- 
tion of transverse wave number q and the initial energy e 
for weak driving U <C Hoj at photon energy ftw = 2meV ~ 
500GHz. The dashed lines separate regions with propagat- 
ing and evanescent modes under the barrier, (c) Cut at zero 
energy revealing the behavior in the evanescent region. The 
dashed line marks the static transmission 2T for graphene 
at energy hw, while the dotted curve is its semiclassical ap- 
proximation. The transverse wave number q is scaled by 
qu, which is defined by q^L — fiw/E^ for graphene, and 
(qujL) 2 = Hoj/El for the 2DEG. 



in the barrier region. Whether modes are evanescent or 
propagating furthermore depends on their transverse mo- 
mentum q. Their response to driving is encoded in the 
function AT/p, which represents the differential response 
for a given mode with transverse momentum q and total 
energy e [see Eq. ((16])]. 

Figure [4] depicts AT/p for graphene and the 2DEG. 
In the 2DEG case, at each energy only a discrete set of 
resonant modes contributes to the pump effect. These 
resonances correspond, up to a shift of ±Huj in energy, to 
quasibound levels in the static system, that form because 
of the velocity mismatch at the interface with the metal- 
lic leads, which results in strong confinement. The con- 
tact resistance increases with increasing Fermi momen- 
tum and the resonances become increasingly nar- 
row, resulting in a suppressed response of the 2DEG. By 
contrast, in graphene a broad range of modes contributes 
at all energies, and the response is particularly strong in 
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FIG. 5. Evanescent mode pump mechanism. Evanescent 
modes penetrating the barrier from the left absorb or emit 
photons such that they become propagating. The correspond- 
ing modes from the right lead decay before reaching the region 
with the ac gating. Filled-in arrows indicate strong popula- 
tion; no n- filled arrows mark negligible scattering channels. 



a diamond-shaped region within the bipolar regime (red 
square in Fig.j4|, in agreement with our earlier obser- 
vation in Fig. |3f Graphene's response does not exhibit 
sharp resonances due to the fact that, in the static case, 
carrier chirality prohibits their confinement, so that carri- 
ers in the graphene pump remain strongly coupled to the 
leads, even when they are driven far out of equilibrium. 
Instead, graphene's response in Fig. [4ja) has features as- 
sociated with the cone e = Hvpq (dashed lines) and its 
replicas shifted by multiples of Huj. For the weak driving 
considered here, only the two replicas at e = Hvpq ± Huj 
are visible. The first cone divides the modes entering 
the scattering region into two categories, propagating 
and evanescent: any incoming carrier in mode q will 
become evanescent under the barrier if its energy ful- 
fills |e| < Hvjrq : and will remain propagating otherwise. 
The other two cones determine whether carriers popu- 
late evanescent or propagating modes after absorption or 
emission of a single photon. Analogously, in the 2DEG 
the threshold between propagating and evanescent modes 
is found at e = ft 2 g 2 /2m*, and under photon emission or 
absorption shifts to e = h 2 q 2 /2m* ± Huj. 

In both systems, it can now be seen that the main 
mechanism of charge transfer is established by a process 
whereby an evanescent mode coming from the left lead 
may get transmitted to the right by absorbing or emitting 
a photon, as long as in this process it jumps to an existing 
propagating mode that may travel into the right lead. 
An evanescent coming from the right, in contrast, first 
encounters the static region not covered by the gate, and 
so gets reflected with a high probability (see Fig. 5j. This 
spatial asymmetry rectifies the carrier flux such that a 
net transport from left to right emerges. 

In the 2DEG, when photon absorption is at resonance 
with a quasibound state, this process via evanescent 
modes yields a maximal value AT/p ~ 1 5 but because 
the resonances are narrow the integrated current remains 
small. In graphene, on the other hand, all evanescent 
modes with qL < hcj/E^ exhibit a response AT/p < 2, 
close to twice the optimal value of resonant evanescent 



modes in a 2DEG, but without the need to satisfy any 
resonance condition. The factor two comes from the two 
possible transitions to a propagating mode, by absorp- 
tion but also by emission of a photon, i.e., it is directly 
related to bipolarity. Resonant conditions are not re- 
quired because Klein tunneling in the propagating modes 
and macroscopic tunneling in evanescent modes keep the 
contacts always open; this is directly related to chirality. 
In combination, these two graphene-specific effects max- 
imize the differential pump response around the Dirac 
point, and thus determine the robust characteristic fea- 
tures of the total pump current in Fig. [3j 



V. CONCLUSIONS 

With this work, we have put forward a pump mech- 
anism that makes use of the evanescent modes under a 
potential barrier in graphene. If an ac gate voltage acts 
upon the, say, left half of the barrier, electrons penetrat- 
ing the barrier from the left lead are excited into prop- 
agating modes and, thus, will be scattered to the right 
lead. By contrast, electrons from the opposite side will 
reach the driving region only with exponentially small 
probability. This breaking of spatio-temporal symme- 
tries induces surprisingly large non-adiabatic pump cur- 
rents in the range of several nA, significantly larger than 
what has been observed so far with semiconductors. The 
main reason for this efficiency is that a whole range of 
evanescent modes is excited, while in the corresponding 
setup with a 2DEG, only isolated resonances contribute. 
Despite the large pump currents, the main effect stems 
from single-photon absorption, which allows one to ob- 
tain analytical results. Moreover, owing to current con- 
servation, ratchets which are built by several pumps in 
series, exhibit identical operation characteristics (repeti- 
tion of the element is therefore not desirable for practical 
implementations) . 

A further important observation is that the resulting 
pump current changes smoothly, symmetrically, and al- 
most monotonically with the distance of the Fermi en- 
ergy to the Dirac point. This behaviour is a direct con- 
sequence of bipolarity and, thus, particular to graphene 
pumps. As an application, it allows steering the dc cur- 
rent into a direction of choice by simply shifting the bar- 
rier height via a local dc gate voltage across the Dirac 
point. This effect can also be used to detect a static elec- 
tric field by measuring the current in response to a small 
oscillating probe electric field, or detecting the amplitude 
of such an oscillating field when the static field is fixed 
to a moderately large value, beyond which the response 
flattens out. 
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